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Abstract. The Large Hadron Collider (LHC) experiments have revealed that the 
predictions of the color glass condensate (CGC) tend to underestimate the multiplicity 
at mid-rapidity. We develop and estimate a full second-order viscous hydrodynamic 
model for the longitudinal expansion to find that the CGC rapidity distributions are 
visibly deformed during the hydrodynamic stage due to the interplay between the 
entropy production and the entropy flux to forward rapidity. The results indicate 
the importance of viscous hydrodynamic evolution with non-boost invariant flow for 
understanding the CGC in terms of the heavy ion collisions. 



1. Introduction 

The heavy ion program at Large Hadron Collider (LHC) opened up a new era in the 
physics of the quark-gluon plasma (QGP) pQ at higher energies. One of the most unique 
properties of the hot matter is the near-perfect fluidity, which was first discovered 
in the Au-Au collisions at Relativistic Heavy Ion Collider (RHIC) [2]. A standard 
modeling of the heavy ion collision at RHIC consists of several stages; nucleus-nucleus 
collision, early thermalization, hydrodynamic evolution, freezeout and hadronic cascade. 
The relativistic hydrodynamic model describes the intermediate stage (r ~ 1-10 fm/c) 
where the system is locally equilibrated. On the other hand, the pre-collision state is 
considered to be described by the color glass condensate (CGC) where the medium is 
interpreted as saturated gluons [3]. The CGC itself is considered to be successful in 
reproducing the rapidity distributions and the multiplicities observed at RHIC. The 
latest Pb-Pb collisions at LHC, however, revealed that most of the CGC predictions 
underestimated the multiplicity at mid-rapidity [4J. This could be due to the fact 
that secondary interactions are missing in the comparison. In this study we estimate 
the viscous hydrodynamic evolution of the CGC rapidity distributions at RHIC and 
LHC with non-boost invariant flow for the first time [5]. We solve the full second 
order constitutive equations [6] with both shear and bulk viscosity to investigate the 
hydrodynamic deformation of the initial distributions. 
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2. Viscous Hydrodynamic Model 

The viscous hydrodynamic equations consist of the conservation laws and the 
constitutive equations. In the limit of vanishing chemical potentials we consider, the 
former is the energy- momentum conservation d^T^ u = 0. We employ the full second- 
order extended Israel- Stewart theory P, [7j for the latter. We choose the Landau frame 
where the flow w M is identified as the eigenvector of T^ v . Then the non-equilibrium 
dynamics in the hydrodynamic system is described by the bulk pressure II and the 
shear stress tensor tc^. The constitutive equations for the dissipative currents read 



where D = u^d^ and V M = <9 M — u^D are the time- and the space-like derivatives. T 
is the temperature and rj, (, r and \ are the transport coefficients. The angle brackets 
on indices (...) denote the projection of the traceless symmetric components. We solve 
the hydrodynamic equations in the (l+l)-dimensional relativistic coordinates to discuss 
non-boost invariant evolution of the medium, although the effects of the transverse flow 
are neglected. This is in good contrast to the (2+l)-dimensional viscous hydrodynamic 
models which are used to discuss transverse dynamics but with boost-invariant flow. 

One has to introduce the equation of state and the transport coefficients as inputs 
to perform hydrodynamic calculations. Here we employ the equation of state from 
the latest lattice QCD calculations [5]. Since the transport coefficients are difficult to 
calculate in the first principle method, the shear viscosity is set to r\ — s/4ir [5] and 
the bulk viscosity ( = |(| — c 2 s )n by extending the method in Ref. [10| where c s is the 
sound velocity. The second order transport coefficients r and x are calculated with the 
ratios of the first and the second order ones given in kinetic theory. Note that they are 
employed for the purpose of demonstration to see qualitative responses of the systems. 

The initial energy distributions are calculated from the Monte-Carlo version [TT|[T2"] 
of Kharzeev-Levin-Nardi (MC-KLN) model [T3| HI] . The saturation scale Q s for a 
nucleus A at a transverse coordinate x± is set as 



where x is the momentum fraction of incident partons and Ta(x±) the thickness function. 
The parameter A is related with the rapidity dependence of the distribution; dN^/dy 
gets steeper with increasing A. Here it is fixed as A = 0.28 [15J. The 5% most central 
events are used for the construction of the smooth initial conditions. The initial energy 
distribution eo(T ,?7 s ) is obtained from the average CGC transverse energy distribution 




(1) 






(3) 
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per unit area {1/ S^^dEx / dy by identifying momentum rapidity y with space-time 
rapidity rj s . Here the initial time is set to r = 1 fm/c. The initial distributions for the 
dissipative currents are not well-known. Here we choose U.(tq, r] s ) = and ^"(to, tj s ) = 
to employ the same initial energy-momentum tensor for ideal and viscous hydrodynamic 
cases. The initial flow rapidity is set to Yf(r , r) s ) — 7] s where Yf — | In rgjj^g • 

We investigate the rapidity distribution through the entropy distribution per flow 
rapidity dS/dYf, where the two quantities are related as dN^ dro /dy ~ (2/3) x (1/3.6) x 
dS/dYf [16J. This is because the ratio of the number density to the entropy density is 
temperature dependent in the limit of relativistic massless gas, and flow rapidity in a 
fluid element on average can be identified with momentum rapidity. Note that this is 
a pure hydrodynamic quantity which does not involve complicated model assumptions 
such as freezeout and thus one can compare the initial and the final distributions. 



3. Results 




Figure 1. The initial rapidity distributions of the CGC (solid lines) modified by ideal 
(dotted lines) and viscous (dashed lines) hydrodynamic evolution. (Left) The RHIC 
case with the final time r/ = 30 fm/c. (Right) The LHC case with 77 = 50 fm/c. 



The CGC rapidity distributions with and without hydrodynamic evolutions for 
Au-Au collisions at ^snn = 200 GeV and Pb-Pb collisions at ^snn = 2.76 TeV are 
shown in Fig. [TJ They correspond to the RHIC and the LHC settings. The distributions 
are always flattened in ideal hydrodynamic systems due to the entropy flux to forward 
rapidity driven by the pressure gradient in rapidity direction. On the other hand, 
the distributions in viscous hydrodynamic systems are enhanced due to the entropy 
production. The entropy production could be larger because the viscous coefficients here 
are close to the conjectured lower boundaries. The final times are chosen so that the 
temperatures at mid-rapidity are sufficiently near the QCD pseudo-critical temperature. 
For the current parameter settings, the equal-time surface is not so different than the 
isothermal one because the distributions do not change much after r ~ 20 fm/c. 



Viscous Hydrodynamic Evolution with Non-Boost Invariant Flow for the Color Glass Condensate^ 

The hydrodynamic modifications differ in magnitude for the RHIC and the LHC 
cases. If, as indicated in the current parameter settings, the CGC rapidity distribution 
is flattened at RHIC, the actual initial distribution needs to be steeper, i.e., the true A 
has to be larger than the current value. Since the hydrodynamic effects are accidentally 
cancelled at LHC, it translates into larger multiplicity at LHC. This could be one of the 
candidates to explain the underprediction of the multiplicity at LHC by the CGC. 

4. Conclusions 

We developed the theoretical and numerical scheme of the full second order viscous 
hydrodynamic model for the longitudinal expansion. The CGC initial rapidity 
distributions are visibly modified during the hydrodynamic stage due to the interplay 
between the two factors: (i) entropy production from non-equilibrium processes and (ii) 
entropy flux to forward rapidity caused by non-boost invariance. The result implies 
that readjustment of the CGC parameters is necessary, which could be one of the key 
factors for explaining the gap between the CGC predictions and the LHC data [4] . The 
results indicate that non-boost invariant hydrodynamic evolution together with viscosity 
is indispensable for developing an integrated picture of the heavy ion collisions. Detailed 
analyses on the parameter dependences will also be explored elsewhere [17]. 
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